Draft version February 1, 2008 

Preprint typeset using I^T^X style cmulateapj v. 10/09/06 



o 
o 

(N 

a 
>—> 

Oh; 

6 : 



> 

O 
o 
p 

o 
i> 
o 



X 



THE SPATIAL STRUCTURE OF AN ACCRETION DISK 

Shawn Poindexter 1 , Nicholas Morgan 1 , Christopher S. Kochanek 1 

Draft version February 1, 2008 

ABSTRACT 

Based on the microlensing variability of the two- image gravitational lens HE 1104-1805 observed 
between 0.4 and 8 /jm, we have measured the size and wavelength-dependent structure of the quasar 
accretion disk. Modeled as a power law in temperature, T oc i? _/3 , we measure a B-band (0.13 /zm in 
the rest frame) half-light radius of -R1/2.B = 6.7l|; 2 X 10 15 cm (68% CL) and a logarithmic slope of 
= 0.6112$ (68% CL) for our standard model with a logarithmic prior on the disk size. Both the scale 
and the slope are consistent with simple thin disk models where (3 = 3/4 and R1/2.B = 5-9 X 10 15 cm 
for a Shakura-Sunyaev disk radiating at the Eddington limit with 10% efficiency. The observed fluxes 



favor a slightly shallower slope, [3 = 0.551q'q2, and a significantly smaller size for /3 = 3/4. 

Subject headings: accretion, accretion disks — gravitational lensing — quasars: individual (HE 1104- 
1805) 



1. INTRODUCTION 

A simple theoretical prediction for thermally radi- 
ating thin accretion disks well outside the inner disk 
edge is that the temperature diminishes with radius as 
T cx i?" 3 / 4 (|Shakura fc Sunvaevl [i~973l ) . This implies a 
characteristic size at wavelength A of R\ oc A 4 / 3 , where 
R\ is the radius at which kT = hc/X. Ne edless to say , 
it is unlikely that disks are this simple (e.g. lBlaesll2004r i. 
but measurement of the size-wavelength scaling would 
be a fundamental test for any disk theory. While the 
angular sizes of quasar accretion disks are far too small 
to be resolved by direct observation, gravitational lenses 
can serve as natural telescopes to probe accretion disk 
structure on these scales. 

Here we make the first such measurement of this size- 
wavelength scaling. Our approach, gravitational mi- 
crolensing of a quasar, will be unfamiliar to the AGN 
community, but it is the only method with the necessary 
spatial resolution a vailable to us for the foreseeable future 
(see the review by Wambsganss 2006]) . The stars in the 
lens galaxy near the image of a multiply-imaged quasar 
generate complex caustic networks with a characteristic 
scale called the Einstein radius Re set by the mean stel- 
lar mass (M). For the lens system we consider here, 
HE 1104-1805, (R E ) = 3.6 x 10 16 ((M) /HMq) 1 / 2 cm. 
Because the magnification diverges on the caustic curves 
of the pattern and the source is moving relative to the 
lens and observer, microlensing allows us to study the 
spatial structure of anything smaller than (Re)- Us- 
ing microlensing to probe scales smaller than Re has 
been successfully appli ed to resolve stars i n Galactic mi- 
crolensing events (e.g. lAlbrow et ani2001h . For quasars, 
this has b een considered analytically or with simula- 
tions (e.g. lAgol fc Kroliklfi99l iGoicoechea et al.l l200l 
iGrieger. Kavser fc Schrammlll991f L but the data and al- 
gorithms needed to implemen t the programs have only 
become available recently (see lKochanek|[2004l) . 

HE 1104-1805 is a doubly imaged radio-quiet quasar 
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at z„ = 2.319 with a separation of 3'.'15 (Wisotzki et al. 
1993). Th e lens at zi = 0.729 was discovered in the 
near-I R bvlCourbin. Lidman fc Magainl (119981 ) and with 
HST (|Remv et al.lll998ULehar et al.ll2000D . Here we an- 
alyze 13 years of photometric data in 1 1 bands from the 
mid-IR to B-band using the methods of lKochaneld (|2004l ) 
to measure the wavelength-dependent structure of this 
quasar modeled as a power law R x oc X 1 ^. We assume 
a flat ACDM cosmological model with = 0.3 and 
Hq = 70 km s _1 Mpc -1 and report the disk sizes assum- 
ing a mean inclination of cosi = 1/2. In $2] we describe 
the data set and our methods, and SJ3] presents our mea- 
surement results and our conclusions. 



2. DATA AND METHODS 
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We included observations of HE 1104-1805 in 
bands: B, V, R, I, J, H, K, 3.6 iim, 4.5 fim, 5.8 /zm, 
and 8.0 fj,xn. These included our own SMARTS 
optic al /near-IR, SOAR near- IR, HST and Spitzer IRAC 
dat a (IPoindexter et al.l|2007t ). R-band monitoring data 
by lOfek fc Maozl (120031). the V-ba rid monitoring data 
from ISchechter et all ( 20031 ) and IWvrzvkowski et all 
(l2003f). and earlier data from iRemv et all ( 1998[) . 
Gil-Merino. Wisotzki fc Wambsganss! (|2002ft . 
Courbin et all (|1998l) . and iLehar etail (|2000l) . Where 
possible we corrected the light curves for the 152 
day time delay betwe en the images we measured in 
iPoindexter et all ((2007). Where we could not, we broad- 
ened the photometric uncertainties by 0.07 mag so that 
the flux ratio uncertainties would be larger by the 0.1 
mag shifts we found between time-delay corrected and 
uncorrected flux ratios. The light curve is plotted along 
with one of the light curve models f rom our ana l ysis i n 
Figure Q] As pointed out in IPoindexter et al.1 (|2007h . 
image A has slowly switched from being bluer than 
image B to being redder in the optical/near- infrared 
(see Figured]), while the mid-infrared fl ux agrees with 
the fl ux ratio of the broad emission lines (Wisotzki et al. 
fl993h . 

Determining the structure of the disk as a function of 
wavelength from such data is relatively straight forward. 
The divergences on the caustics of the microlensing mag- 
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Fig. 1.— The HE 1104-1805 multi-band light curves. The curves 
show a model for the BVRIJHK bands and the four IRAC channels 
(3.6 to 8.0 fim). For clarity we split the data into three panels and 
show the B, V, R, and J-band data points in 200 day averages. 

nification patterns are only renormalized by the finite size 
of the source quasar because the observed magnification 
is the convolution of the pattern with the source struc- 
ture. Thus, larger emission regions will show smaller 
variability amplitudes than smaller emission regions be- 
cause they smooth the patterns more. As we go from the 
K-band to the B-band, the radius of a standard thin disk 
with R\ oc A 4 / 3 , changes by a factor of 8.6, corresponding 
to a change in the disk area of almost two orders of mag- 
nitude. We see in Fig. [l] that the bluer wavelengths show 
larger amplitudes than the red wavelengths, so we im- 
mediately know that the blue emission regions are more 
compact than the red. 

Our analysis use s the Bayesian Monte Carlo method 
of iKochanekl ((2004) to analyze the data. In essence, we 
randomly draw large numbers of trial light curves from a 
range of physical models, fit them to the data and then 
use Bayesian statistics to derive probability distributions 
for the disk structure. We need, however, to discuss the 
physical variables used in the models over which we aver- 
age as well as our model for the structure of the accretion 
disk. 

We use the lens model sequence from lPoindexter et al.l 
(|2007l ). which consists of a de Vaucouleurs model 
matched to the HST observations embedded in an NFW 
halo, to set the shear 7, convergence k, and stellar frac- 
tion k*/k for the microlensing magnification patterns. 
The models were constrained to match the flux ratios of 
the mid-IR IRAC bands. We used a mass function of 
dN/dM oc A/ -1 - 3 with M max /M min = 50 whose struc- 
ture is bro adly consisten t with the Galactic disk mass 
function of iGouldl (|2000f ). The lens models were pa- 
rameterized by /m/Lj the fraction of a constant mass- 
to-light ratio (M/L) represented by the visible galaxy. 
For each of ten models, Jm/l — 0.1, 0.2, . . . , 1.0, we pro- 
duce 2 magnification patterns with an outer dimension 
of 10.Re and 8192 2 pixels to achieve a pixel scale of 
4.4 x 10 13 ((M) //lAfg) 1 / 2 cm/pixel that is smaller than 
the gravitational radius r g — GM/c 2 = 3.5 x 10 14 cm 
expected for HE 1104-1805 (see ij3]). We experimented 
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Fig. 2. — Probability distributions for the B-band accretion disk 
half-light radius. The bold solid (dashed) curve is the distribution 
based on a logarithmic (linear) prior on the disk size. The filled 
(open) square near the top is the median of the distribution with a 
logarithmic (linear) prior along with the 68% confidence error bars. 
The lower curves show the contribution from the more significant 
j3 trials, labeled by their value of f3. The vertical lines show the 
B-band thin disk size predicted from the I-band flux (Eqn. |3j 
and the size expected from standard thin disk theory (Eqn. [2} for 
Eddington limited (L/Le = 1) accretion with r\ = 10% efficiency. 

extensively with magnification patterns of different sizes 
and scales to ensure that our choice of pixel scale and 
outer dimension di d not affect our r esults. We used the 
velocity model from lKochane k (2004) , where the adopted 
values are 73 km/s for our velocity projected onto the 
lens plane, 308 km/s for the velocity dispersion of the 
lens estimated from fitting an isothermal lens model and 
RMS peculiar velocities of 135 and 71 km/s for the lens 
and source respectively. We assume a uniform prior on 
the mean masses of the stars of 0.1 Mq < (M) < 1.0 M , 
but this has only limited effects on the estimates of the 
disk size (see Kochanek 2004) . Our estimates of the disk 
properties average over all these parameters. 

We modeled the disk as a face-on thermally radiatin g 
disk with temperature T oc Rr 13 (e.g. lCollier et al.1l!998h . 
corresponding to a surface brightness profile of 



/„ oc v 3, exp (R/R) 
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(1) 



where the size scale is R\ = i?e(A/ 'Xb) 1 ^ and Rb is 
the disk size at the observed B-band (1310A in the rest 
frame) and (3 = 3/4 for standard thin disk theory. In 
many cases we report the half-light radius R\/2,\{P) — 
C(f3)R x (/3) where for = 3/4, C = 2.44. While we 
do not include an inner disk edge of Ri n ~ 2r g to 6r g , 
the lack of this central hole has little effect on our re- 
sults unless /3 — > 2 or Rb — Rin- We tried disk models 
with profile exponents of j3 = 1/4, 3/8, •■ • ,7/4, making 
4 x 10 6 trial light curves for each value of /3. In our final 
analysis we use the 3.3 x 10 6 light curves that passed a 
threshold of x 2 /Ndof < 3 for Ndof — 438 degrees of free- 
dom. We used both a logarithmic prior, P(R\) oc l/R\, 
and a linear prior, P(R\) oc constant, on the disk size. 
Generally, a logarithmic prior is preferred for scale free 
variables like R\ , but for this problem a linear prior may 
be more appropriate for small source sizes because we are 
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Fig. 3. — B-band half-light disk radius i?i/2,B versus {3 where 

T oc R-t 3 . The solid (dotted) contours represent the 68%, 90%, and 
95% confidence levels of the microlensing measurements assuming 
a logarithmic (linear) prior on disk size. The small solid contours 
is the B-band flux estimate considering a T oc R-P model fit to 
all the bands. The dashed line shows the size estimated from just 
the B-band flux measurement without the constraint on (3 from the 
other bands. The large solid square at j3 = 3/4 shows the radius 
predicted by thin disk theory (Eqn. [2j. 

sensitive to the difference between smali sizes oniy dur- 
ing caustic crossings. If, however, we have solutions with 
the caustic crossing sitting in a gap of the light curve, 
we will find that P{D\R\) approaches a constant value 
as R\ goes to zero. This leads to a formal divergence 
in P(R\\D) with a logarithmic prior, suggesting that a 
linear prior may be more appropriate. In practice, Figs. 
[2K3] show the results for both options and we use the re- 
sults for the logarithmic prior in our discussion. For the 
results we give the value at the median of the probability 
distribution and the 68% (1 a) confidence regions. 

3. RESULTS AND DISCUSSION 

We have no difficulty reproducing the observed light 
curves, including the presently observed color reversal. 
We illustrate this in Fig. [TJ where we superpose our best 
fitting light curve model on the data. 

We find a wavelength-size scaling of the accretion disk 
of P = 0.61±g i7 and f3 = 0.89t^ for the logarithmic 

and linear priors respectively (Fig. 0J, where R\ oc X 1 /^ 
and T oc R~P. This is consistent with simple thin 
disk theory (J3 = 3/4). The B-band half light radius 
is R\/2,b = 6.7^3 2 x 10 15 cm with the logarithmic prior 



and -Ri/2,b = 1-3 



+0.9 
-0.6 



10 cm with the linear prior 



(Figs. [2] and [3|). We use the half light radius rather 
than R\ because it has less covariance with the expo- 
nent (3. We can also estimate the sizes for the individ- 
ual bands, as shown in Fig. [4j although these will be 
highly correlated because the model only has the expo- 
nent j3 and one scale length as actual parameters. If we 
fix (3 — 3/4, then we find that B-band size for this case 

is Ri/2,b(P = 3/4) = 9.0t^ x 10 15 cm. 

We can compare to thin disk theory only for the case 
of (3 — 3/4 since there is no simple generalization for 
alternate values of (3. Based on the CIV emission line 
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Fig. 4. — Half-light disk radius R\/2,\ versus wavelength. The 
filled squares (open triangles) are the size estimates from microlens- 
ing (the source flux, Eqn. |3j. Note that the microlensing uncer- 
tainties are highly correlated because the only 2 actual variables 
are i?B an d 0. The short dashed and dotted lines are the best 
fit power laws to these measurements. The long dashed line shoes 
how the normalization of the points would shift if we use j3 = 3/4 
rather than the best fit slope of /3 = 0.55. 

width of HE 1104-1805. [Peng et all (|2006l ) estimated a 
black hole (BH) mass of Mno4 = 2.4 x 1O 9 M , which 
corresponds to a gravitational radius of r g = GM/c 2 = 
3.5 x 10 14 cm. In thin disk theory (|Shakura fc Sunvaevl 
Il973f ). this implies a size scaling of 
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hp is the 



where rj is the radiative efficiency (L = r/Mc 2 
Planck constant, and L/Lr is the fraction of the Edding- 
ton luminosity radiated by the QSO. Thus the expected 
B-band (rest frame 0.13 jiva) size for a disk radiating 
at the Eddington limit (L/Le) with 10% efficiency is 
Rb — 2.4 x 10 15 cm corresponding to a half-light radius 
of R\/2,b — 5.9 x 10 15 cm. This agrees well with our 
measurement (Figs. [21 El and [J). For the (3 = 3/4 model 
the disk scale length is significantly larger than the grav- 
itational radius (i?i/2,B ~ 20r ff ) and corrections for the 
inner edge of the disk will be modest. Particularly if we 
allow for uncertainties in the black hole mass estimate 
(~ 0.3 dex), there is good agreement with the simplest 
possible thin disk model. 

The magnification-corrected flux of the quasar provides 
a second comparison scale under the assumption that the 
disk is thermally radiating. If the quasar has a magnifi- 
cation corrected magnitude of m, then 
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where X b s is the observed wavelength, zpt is the filter 
zero point (normalized to the 1-band), and 



K(J3) 



1 



2.58 



udu [exp(u^) — l] 



(4) 



is the /3-dependent term due to the temperature pro- 
file normalized so that K(J3 = 3/4) = f. Assuming 
the magnifications of images A and B are 11.5 and 4 
respectively (as found in the best fit macro model in 
iPoindexter et al.ll2007h . we calculated the magnification 
corrected source magnitude in each of the ff bands 
to estimate the disk size versus wavelength (Fig. [4]). 
For the H-. I-, and V - band , we used the HST obser- 
vations of iLehar etall (l2000fl. Th e mid-IR magnitudes 
are from Svitzer (|Poindexter et al.|[2007l ). We calibrated 
the SMARTS B/R and J/K data using the Guide Star 
Catalog and 2MASS respectively. 

The sizes estimated from the flux are very well fit 
by a power law (see Fig. [4} with a slope equivalent 
to (3 = 0.55^q'q2 when we assume a magnification un- 
certainty of a factor of 2. While this slope is con- 
sistent with our microlensing results, the size scale of 



R 



i/2,: 



1.8_Q 5 x 10 5 is smaller by a factor of 4 than 
our standard estimate. This discrepancy depends on the 
value of 0, with a flatter temperature profile showing 
less of a difference. While there is little extinction in the 
lens, 2.9 mags of B-band (0.3f /im rest frame) extinction 
in the source would reconcile the flux and microlensing 
size es timates. However, w ith such a n extinction, neither 
SMC dGordon et al.ll200l or AGN ([Gasket! et al.ll200l 
extinction curves can reconcile the two estimates at all 
wavelengths. Nonetheless absorption in the source could 
be a partial explanation. 

The general relationship we find between the mi- 
crolensing, thin dis k and flux sizes seems to be typical. 
IPoolev et al.l ([2007) noted qualitatively that microlens- 
ing sizes tended to be larger than expected from t he op - 
tical flux and thin disk models. iMorgan et al.l (|2007l ) 
showed quantitatively that the microlensing sizes scaled 
as expected with BH mass (i? A oc M 2 / 3 ) and were con- 
sistent with being proportional to the flux sizes, but that 
the absolute scales of the microlensing sizes were slightly 
larger than the thin disk sizes and considerably larger 
than the flux sizes. Our results here suggest that part of 
the solution may be that the effective temperature profile 

is somewhat shallower than T oc R~ 3 / 4 . 

Several local estimates (|CoHier et al.l 119991 : 
iSergeev et al.l 120051 ) have found UV-optical wave- 
length dependent time delays of nearby AGN consistent 
with T oc i? -3 / 4 . They also found the flux discrepancy, 



but phrased the problem as needing to put the systems 
at higher than expected distances (through a low value 
of Ho) in order to reconcile the model disk surface 
brightness with the observed flux. 

Our current disk model is a face-on, thermally ra- 
diating disk without the central temperature depres- 
sion created by the inner edge of the accretion disk 
(Eqn [1} . Omitting the inner edge has little effect be- 
cause at fixed wavelength little flux is radiated there 
and the finite resolution of the magnification pat- 
terns eliminates the formal surface brightness diver- 
gence in Eqn. |T]). It is expected from theoretical 
work that microlensing is primarily sensitive to an ef- 
fective smoothing area, and the resulting size estimate 
is only wee kly sensitive to the true surface brightness 
profile (e.g. iMortonson. Schechter fc Wambsganssl 120051 : 
ICongdon et al.ll2007l ). Nonetheless, a clear next step is to 
begin interpreting the micro lensing results using more re- 
alistic thin disk models (e.g. Hub env et al.ll2001l ) to see if 
adding such details alters the basic picture. The present 
results suggest that part of the solution may be to al- 
ter the temperature profile of the disk, perhaps through 
irradiation of the outer regions by the inner regions. 
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